Smectic-A elastomers with weak director anchoring 



J. M. Adams and Mark Warner 
Cavendish Laboratory, JJ Thomson Avenue, Cambridge CBS OHE, United Kingdom 



Olaf Stenull and T. C. Lubensky 
Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, 



USA 



00 

o 

O 

(N 



«4H 

o 



I 

C 

o 
o 



> 

On 

m 
oo 
p 

o 

00 
O 



X 



Experimentally it is possible to manipulate the director in a (chiral) smectic- yl elastomer using 
an electric field. This suggests that the director is not necessarily locked to the layer normal, as 
described in earlier papers that extended rubber elasticity theory to smectics. Here, we consider the 
case that the director is weakly anchored to the layer normal assuming that there is a free energy 
penalty associated with relative tilt between the two. We use a recently developed weak- anchoring 
generalization of rubber elastic approaches to smectic elastomers and study shearing in the plane 
of the layers, stretching in the plane of the layers, and compression and elongation parallel to the 
layer normal. We calculate, inter alia, the engineering stress and the tilt angle between director 
and layer normal as functions of the applied deformation. For the latter three deformations, our 
results predict the existence of an instability towards the development of shear accompanied by 
smectic-C-like order. 



PACS numbers: 61.30. Vx, 83.80.Va, 81.40.Jj 



I. INTRODUCTION 

Smectic- A (Sm-A) liquid crystal elastomers incorpo- 
rate the anisotropic properties of liquid crystals, and the 
rubber elasticity of polymer networks. The formation of 
a smectic layer structure by the mesogens is the cause 
of their particularly anisotropic elastic properties. These 
elastomers have been synthesized, and their elastic prop- 
erties explored through mechanical testing. Nishikawa 
and Finkelmann found that a class of strongly coupled 
Sm-^4 systems behave like 2D elastomers in the layer 
plane [l[ but that they are extremely stiff when stretched 
parallel to the layer normal. At a threshold of a few 
percent strain along the layer normal, the elastomer be- 
comes mechanically softer, and turns opaque Q because 
of layer rotation. This behavior is reversible when the 
strain is removed. Other weakly coupled smectics also 
have thresholds, but not to layer rotation, and do not 
have the same extreme mechanical anisotropy Q. A 
threshold has also been reported in Sm-^4 materials with 
a shorter correlation length of the smectic layers, but 
here the sample remains transparent after the threshold 
[J]. Typically, smectic elastomers are synthesized in the 
form of films, either with the layer normal in the plane of 
the film P, 0, [H or with the layer normal perpendicular 
to the film plane However, mechanical testing has 
only been performed on the first of these two types. The 
second type would be useful for performing mechanical 
compression tests parallel to the layer normal. 

Experimental study of the electroclinic effect in Sm-A 
elastomers suggests that the layer normal and the di- 
rector can be manipulated with an electric field and indi- 
cates that the two are not necessarily rigidly locked Q , at 
least on the scale of electrical energies, as a Lagrangian 
elasticity theory developed in Ref. Q assumes. Since 
rubber elastic energies are typically larger than those of 
electric fields, one would expect mechanical fields also to 



induce relative rotations. In fact, such relative rotation 
has been observed experimentally in very recent shear 
experiments by Kramer and Finkelmann [8(. 

On the theoretical side, a model of Sm-A elastomers 
has been constructed using non-linear rubber elasticity 
extended to smectics @ that describes the results of 
Nishikawa and Finkelmann well. This model rigidly locks 
the director to the layer normal. As mentioned above, 
a model based on Lagrangian elasticity theory has also 
been developed ml- This model fits well with the data, 
however, unlike [9(, it allows in principle for the rela- 
tive rotation of the director and the layer normal. Trig- 
gered by experiments of Kramer and Finkelmann 0, [l(| 
where in-plane shears were applied to Sm-A elastomers, 
the rubber elastic approach has been extended I 111! very 
recently to the case of soft-anchoring. In Ref. [LP ], we 
used this model specifically to study the shear experi- 
ments of Kramer and Finkelmann. In the present paper, 
we employ this model to study various shear and stretch 
deformations. 

The plan of presentation is the following: First, we 
review the extension of non-linear rubber elasticity the- 
ory to smectics. The fundamental distortions of imposed 
in-plane stretch and in-plane shear are then explored. 
Armed with these modes of deformation, we then ex- 
plore imposed extension and compression along the layer 
normal which are complex but decomposable into the 
fundamental modes. 



II. MODEL FREE ENERGY 

The model free energy reviewed here is generalization 
of the original rubber elastic model of smectics @ , and 
was developed originally in Ref. [ll|. It has contribu- 
tions from the background ncmatic elasticity, and from 
the compression or dilation of the layers. In addition to 
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these two terms, we include here a potential that penal- 
izes the deviation of the director from the layer normal 
El: 



/ /trace H~ /layer ~\~ /tilt- 



(i) 



The nematic component of the free energy density has 
been widely discussed [l3| , and is given by 



/tra 



A • £o ■ A • r 



(2) 



where /i is the isotropic state shear modulus, A is the de- 
formation gradient tensor, |o = §_+(r~ l)nono is the step- 
length tensor before the deformation has been applied, 
£ _1 = §_ + (l/ r — l)nn is the inverse of the step-length 
tensor after the deformation has been applied, and §_ is 

the unit tensor. The step length-tensor is proportional 
to the second moment of the Gaussian distribution of 
anisotropic chains making up the rubbery network. The 
anisotropy of the chains is parameterized by r. The smec- 
tic layers embedded in the elastic matrix give rise to the 
following contribution to the free energy 
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where B is the layer spacing modulus (typically larger 
than in liquid smectics), d is the current layer spacing, 
do is the natural layer spacing and is the angle between 
the layer normal and the director. For smectics where the 
layers are strongly coupled to the matrix, changes in layer 
spacing can be derived by analyzing how the embedded 
layer normal deforms in step with the elastic matrix. The 
layer normal k and layer spacing are given by 0] 
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where ko is the initial layer normal. The deformation of 
the layer normal outlined above affinely follows that of 
the rubber matrix because the energetic penalty of not 
doing so scales with the system size in the microscopic 
model described in Q. The cosG term reflects the pro- 
jection of the rods making up the layer spacing contract- 
ing as the director tilts. Note that this term does not 
include the finite thickness of the rods, since as O — + it/ 2 
then d — ► to avoid a diverging energy penalty, and so 
the layer thickness tends to zero. However, physically 
there must be a transition to a constant layer thickness 
as the angle increases, i.e., there must be forces pre- 
venting — > 7r/2 independent of the layer-compression 
term. 

The contribution / t nt = /tilt (sin 0) to the free energy 
density penalizes the deviation of the director from the 
layer normal. To ensure a finite layer thickness, see the 
discussion above, the general form of the tilt contribution 



will be such that /tut (sin 0) — > oo as — > ir/2. When 
expanded to leading order in sin©, it reduces to 



/tilt = h a t sin 2 ■ 



(5) 



where at is a coefficient that vanishes as the A-C tran- 
sition is approached. For simplicity, we will work in the 
following with the simple phenomenological form ([5]). 

Typical values of the constants are B ~ 10 7 Pa, at ~ 
10 5 - 10 6 Pa, and /x - 10 5 - 10 6 Pa. In determining the 
main features of the material properties only the ratios of 
these values are important. Consequently we will denote 



b = B I ' \i and c = at/fx. 



(6) 



Note that the limit c — > oo locks the director to the layer 
normal. On the other hand, when c is small as it would 
be near the AC transition, deviations of from zero 
can be large. Dominating these energy scales is that for 
volume change. The bulk modulus is, as for all rubbers, 
of the order of 10 9 Pa which means that deformations are 
at constant volume, that is dot [A] = 1 rigidly. 

Note that a semi-soft term could also be included of 
the form 



/semi o QJ/xTr 
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(7) 



The value of a can be estimated from either the soft 
plateau in elastomers, or from the electroclinic effect, and 
is typically found to be a < 0.1. It turns out that this 
term does not affect the qualitative features of the me- 
chanic response to the deformations that we consider. 
We will discuss this issue briefly and excmplarily in our 
analysis of sz-shear, but we will not scrutinize the effects 
of the semi-soft term in detail. 



III. EXAMPLE GEOMETRIES 

Three deformation geometries will be considered here, 
and are shown in Fig. [TJ together with the initial director 
and layer normal which are assumed to be along z. 
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FIG. 1: The imposed shear, in-plane elongation and out of 
plane elongation (and compression) geometries that will be 
considered in this paper. 

The model described above has a complicated form 
because incomprcssibility forces the appearance of a co- 
factor of the deformation gradient A- It is possible to 
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eliminate the cofactor dependence by recognizing that 
the director and layer normal move in a 2D subspace, 
which we take to be the xz plane (where z is parallel 
to the initial layer normal). Then all the actual calcu- 
lations can be made using a 2D representation of this 
model outlined in appendix,, 



A. Imposed A xz 

First we examine the imposed shear deformation as 
depicted in Fig. [TJ Here the following will be used 
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n = (0,0,1) 
n = (sin#, 0, cos#). 



(8) 



To ensure incompressible response, dct[AJ = 1, one takes 

A yy = 1/(A XX A ZZ ). For this choice of A, the new layer 

normal is A _T ■ ko = ko/A zz . Layers are unrotated, but 

generally dilated by A zz . The deviation of n from ko, 
that is 0, can thus be identified with the usual 9, the 
rotation of the director. We use 9 until later when there 
is layer rotation and the distinction must be drawn. On 
substituting these expressions into the free energy density 
we obtain the following expression 

2 
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(A zz cos 9 + A xz sin 9) 2 + r(A zz sin 9 — A xz cos ( 



For this A, the additional semi-soft term in the free energy 

is /semi = ^M-^ sm2 H docs not change the qualita- 
tive behavior of the elastomer, so will not be considered 
here. 

This free energy should now be minimized with respect 
to 9, A xx and A zz . It is straight forward to find the min- 
imal value of A xx : 

(10) 
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We now assume that B is much larger than a t and /i so 
that we can approximate A zz ~ cos#. Such an identifica- 
tion is forced by the first term of Eq. © when its coef- 
ficient is large. The leading order response for A xz <C 1 
is then (recalling that c is the reduced angular modulus 
a t //i): 



(r - l)rA x 
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(11) 
(12) 

(13) 



where a is the nominal or engineering stress df/d\ xz . 
Asymptotic analysis can also be performed for large de- 
formations, but it should be remarked that whilst the 
assumptions of very large strains seem unrealistic, the 
results obtained from asymptotic methods are usually ap- 
plicable to a much larger region than anticipated. Here 
we assume that A xz 3> 1 , and A zz <C 1 , and again we have 
to leading order: 
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(14) 
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Shear across an unmoving director has a modulus r/x, as 
for instance inspection of Eq. (|13p in the c 3> 1 limit 
reveals. Here in this unphysical limit, the modulus has 
dropped to fj., indicative of minimal chain extension in 
the gradient direction of the shear. Thus the form of the 
stress indicates that the director rotation is very large, 
evidently ~ n/2. As discussed above, this extreme direc- 
tor rotation is an artifact of the reduction of /tut, which 
in its full form has to suppress the approach of to ir/2 
on physical grounds. 

To illustrate the intermediate features of the model, 
Fig. [5] shows the numerical solution to the minimization 
problem, and Fig. [3] shows an illustration of the defor- 
mation of the elastomer. Since the xx and zz response is 
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FIG. 2: The nominal shear stress, the deformation compo- 
nents A zz and A xx , and the rotation angle 9 of the director 
for an imposed A xz deformation. The dashed (red) line has 
(b, c, r) = (60, 100, 2), the solid (green) line has (60, 1, 2) and 
the dotted (blue) line has (1, 1,2). 

independent of the sign of A xz , on symmetry grounds for 
small imposed shear, A xx — 1 ~ A xz and A zz ~ — A 2 Z . Ro- 
tation does sense the sign and hence 9 ~ A xz . For small 
values of c, 9 saturates at ir/2 for large deformation. By 
contrast, for large values of c the director rotation is sup- 
pressed, in agreement with Q. 

In the shear experiments of Kramer and Finkcl- 
mann [lCj| the applied deformations are similar to 
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FIG. 3: An illustration of the rotation of the director, and 
the sympathetic shears for an imposed A xz deformation, with 
(b, c, r) = (60, 1, 2) as on the solid (green) line in Fig. [2] 



the one that we just discussed. There is, however, the 
difference that these experiments use setups (tilter and 
slider) that pre-set A zz as shear proceeds, and hence free 
relaxation of A zz as above is not possible. A detailed 
analysis of these experiments is given in (Tl| . 



B. Imposed A xx 

The deformation tensor is still that of Eq. ([8|), and 
the free energy takes the same form Eq. © as in the 
previous section. The relaxation behavior of the system 
is sketched in Fig. [4j which illustrates several interesting 
features including a threshold at which rotation of the 
director starts, and a non-monotonic stress-strain curve. 
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FIG. 4: The nominal shear stress, rotation angle, and defor- 
mation tensor components for an imposed A xx deformation. 
The key is the same as in Fig. [2] 
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FIG. 5: An illustration of the rotation of the director, and 
the sympathetic shears for an imposed A xx deformation, cor- 
responding to the solid (green) line in Fig. [4] 

Analytically it is possible to obtain an expression for 



Minimization of the free energy in Eq. ^ 
to A xz results in the following expression 



Axz — 



(r — 1)A ZZ cos 8 sin ( 



with respect 



(16) 



We again consider the case where B ^> at, fi, and so A zz ~ 
cos 8. This leaves only the variable A zz to be minimized 
over in the problem. For large B the threshold occurs 
when there is a minimum at A zz ~ 1 , which results in the 
following condition on A c , the critical value of A xx , 



1) - X 2 c (cr - 1) - r = 0. 



Consequently, the threshold is given by 



A = 



(cr - 1) + y/A(r - l)r + (cr - l) 5 
2F 7 !) 



(17) 



(18) 



This threshold is unphysically large for c> 1 and so it 
would be inaccessible to mechanical experiments in this 
limit. But when r = 2 and c = 1, we have A c = V2, in 
agreement with the numerical results presented in Fig. [4] 
and possibly within the range of experiments. 

The behavior of the system for both small and large 
values of the deformation can again be analyzed. Be- 
fore the threshold, we have 8 = A xz = 0, and A zz = 1. 
Consequently 
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which is the 2D rubber clastic response seen experimen- 
tally in the case where the director and layer normal 
are rigidly anchored, and described theoretically in this 
framework in Q . After the threshold we have for A xx ^> 1 
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This limit is again extreme and unphysical, but useful for 
understanding trends. The result is the same as that for 
a ncmatic rubber that is being stretched perpendicular to 
its director, the director subsequently rotating to be par- 
allel to the applied elongation, the x-direction. The 1/r 
reduction in the effective modulus arises because there is 
an effective xx-elongation of ^fr on director rotation and 
the effective extension with respect to this state is only 
A xx /Vr. 



The induced shear at the threshold A x 



A r increases 



the threshold value of A xx at which the instability starts. infinitely quickly with stretch, as is the case in theory 
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and experiment [14[ for the response at the threshold 
when nematic elastomers are stretched perpendicular to 
their initial director. A symmetry argument shows that 
this must be the case: the instability is insensitive to the 
signs of the angle 9 and the shear A xz , and the stretch 
A xx does not distinguishing among these signs. Thus one 
must have A^ z ~ A xx — A c and 9 2 ~ A xx — A c (or higher- 
even powers of A xz and 9). On taking roots, one has 

±VA XX — A c , that is, singular 



A r and 9 



A xz ~ iv A xx - 
growth at A c . 

The results here clearly show a non-monotonic stress- 
strain relation that is not seen, for instance, in theory 
or experiment for Sm-A elastomers at their instability 
under a A zz extension along their layer normal (see sub- 
section UTTC] below). In the zz extension case, the layer 
normal rotates away from the extension direction to al- 
low in-plane stretch and shear. But in rotating away, 
the layer normal takes the nematic director with it and 
thus towards the contraction diagonal associated with 
the shear that is growing in a singular fashion. For pro- 
late chains, this compression along the director of the 
naturally long dimension of the chain distribution costs 
additional energy This is still a lower energy path than 
that associated with layer extension which would be suf- 
fered if layers did not rotate. By contrast here in the 
case of in-plane stretch, Fig. [5] shows that when director 
rotation takes place, it is instead towards the direction of 
the extension diagonal of the induced xz-shear and thus 
the naturally longer dimension of the chain distribution 
is more accommodated. This shear grows in a singular 
fashion. The consequent slow down in the growth of the 
clastic free energy is sufficient to make the slope of the 
stress-strain relation negative. 

Just as for the van der Waals gas, this type of non- 
monotonic stress-strain curve is mechanically unstable to 
disproportionation. It can be rectified by taking a mix- 
ture of the rotated and unrotated states, and making a 
Maxwell construction on the stress-strain relation. Ex- 
perimentally, it is expected that there will be a plateau 
in mechanical behavior, that is reversible. This is highly 
unusual clastic behavior for a solid. 




FIG. 6: An illustration of the free energy and stress for 
the non-monotonic stress-strain. The curve illustrated has 
(b,c,r) = (60,1,2). Here A A = 1.34, A B = 1.87, and 
cr p = 0.978^. 



existing strain with the same engineering stress <j p . and 
we determine o~ p and /o such that / just touches / = 
at two points. An illustration of / which highlights its 
non-convexity is shown in Fig. together with an illus- 
tration of the deformation gradient at the two minima. 
The stress is also shown, with a plateau, that is formed 
by mixing the two deformation gradients illustrated. The 
anticipated microstructure may be as illustrated in Fig. [71 
however the effects of surface energy and interfacial re- 
gions may lead to a different structure. The actual region 




To carry out the Maxwell construction, we define / = 
/— (o~ p \ xx + /o), he., we subtract from he free energy the 
common tangent that touches /(A xx ) at the points of co- 



FIG. 7: An illustration of a possible microstructure that 
would result from the applied in-plane deformation. Depend- 
ing on the cost of the interface the disproportionation illus- 
trated here may happen in several regions in the sample. The 
region where the microstructure forms may be opaque as il- 
lustrated. 

of the singular shear and director rotation seen in Fig.[?]is 
thus jumped across and these interesting responses may 
not be observable. 

The value of the threshold provides a useful test for 
the value of the modulus at- As no threshold has been 
observed through in-plane elongation experiments, it is 
thought that at should be comparable to or larger than 
\x. For some samples, these stretching experiments have 
been performed up to a strain of A xx = 2.5, but no devi- 
ation from linearity has been observed [loT ] . This suggest 
that the anchoring is so strong in some elastomers, that 
the sample ruptures before the director and the layer nor- 
mal become unlocked. 



C. Imposed A zz 

It is well-known that in liquid Sm-A and in elas- 
tomeric Sm-A Q there is an instability when stretching 
parallel to the layer normal. For small deformations the 
elastomer simply elongates, however for large deforma- 
tions above a threshold, layers start to rotate, and the 
sample undergoes effective in-plane shear, as first pre- 
dicted for elastomers within continuum elasticity [l7| and 
previously analyzed in the framework of rubber elastic- 
ity Q but with n and k rigidly locked (c = oo). When 
c is finite, elongation parallel to the layer normal still 
results in the same instability, as is already known from 
Lagrangian elasticity methods 0]. More unusually, as 
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a consequence of finite c, we predict that compression 
(A zz < 1) can also result in a rotational instability, but 
only of the director with an unrotating layer system. 



1. Elongation (A zz > 1) 

Contrary to the deformations considered thus far, elon- 
gation along the initial layer normal generates layer ro- 
tation, and we are compelled to distinguish between di- 
rector and layer rotation. Thus we take to denote the 
angle of the director with respect to the current layer 
normal, and we introduce £ as the angle by which the 
layers rotate relative to the initial layer normal. The 
deformation gradient tensor 



A xx 











i 







A zx 





A zz 



(25) 



must be assumed since A zx shears are those that produce 
layer rotations, see Eq. This awkward form of the de- 
formation gradient can be decomposed to reveal the true 
deformations by incorporating a rotation matrix. We 
first deform the system by A which we take to be as in 

the previous cases, i.e. as in Eq. ((5|) with stretches and 
instead an xz-shear. The layer normal, given this type 
of deformation, is as yet unrotatcd. We then perform a 
body rotation on both the elastomer and the director, 
which of course leaves the free energy invariant. The ro- 
tation can be chosen to transform the xz shear of A into 



the zx shear of A in (f2"5|) . The overall deformation is 
made up as follows 



A x 



A* 



A = R(0 • 4 ; 4 : 



Ayy 

A zz 



(26) 



where ii(C) represents a rotation through an angle £ 

about y. At this point we can see that all mechanical de- 
formations of Sm-A can be reduced to the three of Fig. [I] 
- others can be reduced to these by suitable rotations. 

Fig. [5] shows the development of the Ay, including A zx 
in response to A zz . The last sketch has the A xz marked in 
as it would appear in the frame of the body. The rotation 
angle ( can be determined by demanding that the A zx 
component is zero. Practically, this is done by putting in 
the explicit form for Ji in the above, evaluating A_ = i2 T A 

and inspecting A zx . It is also clear geometrically from 
the second to last frame of Fig. [8] that to convert A to 

the lab frame A one needs to rotate by the shear angle 
( = tan _1 (A xz /A zz ), resulting in cos£ 



-=4^ . The 

V A -+AL 

connection between the new deformation components Ay 




FIG. 8: An illustration of the rotation of the director, and 
the sympathetic shears for an imposed A zz elongation. The 
A xz shear is shown in the frame of the layer system. £ is the 
shear angle and the rotation angle to go from the lab to the 
layer frames. 



and the old Ay is then 



^A xz + A Z 2 Z 
V A xz + AL 



A 



A 2 



(27) 

(28) 
(29) 



see Q for explicit discussion. 

Recall 9 is the angle between the layer normal and 
the director. Since the body and the system of layers is 
rotated by — £, the director can be expressed as 



n= (sm[e-C],0,cos[8-C]): 



(30) 



so that after rotation by ( in (f2"6"| it becomes 
n = (sin[9], 0, cos[0]). After this transformation is per- 
formed, we obtain the free energy density again given 
by Eq. © with Ay — > Ay and 6 — > 9. However, 
we do not impose A zz , but still impose A zz , whereupon 
A zz = \J A zz — A xz on rearranging ([29)) . The behavior of 
the system can be analyzed numerically and results for 
this analysis are presented in Fig. [9l Note we still see a 
threshold behavior, but now in addition a rotation of the 
director with respect to the layer normal 9. From the 
decomposition discussed above, it is clear the rotation of 
the director with respect to the layer normal arises be- 
cause of the effective xz shear. Consequently we expect 
a non-zero value of 9 on symmetry grounds as observed 

in gnng 

There are several ways that the position of the thresh- 
old can be computed. One method is to expand the free 
energy density for small values of 9. We can then min- 
imize with respect to 9, and substitute in the minimal 
value. It is then possible to minimize over Axx, and sub- 
stitute this back the free energy. Only A xz remains to 
be minimized over. Setting the second derivative of this 
expression with respect to A xz to zero at A xz = then 
gives the equation for the threshold, because it gives the 
point where the level sets change from being convex to 
concave. This results in the following polynomial for the 
critical value A r of A„ 



b\(X c - 1) 2 A^ 
+cr A c [(b — r -+ 



-b(X c 
DA? 



-l)A c 2 (2r-l)-(r- 1)^-1) 
&A 2 - l] = 0. (31) 
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FIG. 9: The nominal shear stress, rotation angle, and defor- 
mation tensor components for an imposed A zz > 1 deforma- 
tion. The key is as in Fig. [2] Note that the dashed (red) 
curve for O lies above zero for A zz exceeding the threshold, 
although this is hard to see because O is very small. 



The terms in c have been collected together here, so that 
it is clear that for c> 1 the same polynomial in square 
brackets as in the locked director case @ is recovered. 
In that case there was only an instability if the reduced 
layer spacing modulus was large enough, b > i — 1. Now if 
the director is no longer locked to the layer normal, finite 
c, then there is naturally an instability for even smaller 
values of b since the lower energy route of rotation is more 
accessible if the director can rotate towards the extension 
diagonal associated with the concomitant shear. 
In the limit B > at,/j for which A c exists, 



2r - cr - 1 + \J{cr + 1) (4r 2 + cr - 4r + 1) 

Ac * 1 + 2b~r • 

(32) 

For r — 2 and b = 60 we have: 



lim A c 

c — >oo 

lim A c 

c^l 



1 + T w 1.033. 
b 

1.028, 



(33) 
(34) 



where the limits have been taken from Eq. (|3T|) . Conse- 
quently the experimental evaluation of this threshold is 
not sensitive to the value of c, and does not discriminate 
between the weak anchoring, and the locked layer normal 
theories. 

For small deformations (before the threshold is 
reached) we have 



(35) 
(36) 





e = o, 




l 


XX 


A zz 


/ = 


iB(Azz-l) 2 + ^ 



(37) 



Before moving on to compression, we find it worth- 
while to compare our results to the available experimen- 
tal findings and other theoretical predictions. The stress- 
deformation curves for B 3> \x are in absolute agreement 
with the experimental curves by Nishikawa and Finkel- 
mann [l[ as well as the theoretical predictions of Refs. @ 
and 0- The results, or at least their interpretation, dif- 
fer, however, as far as O is concerned. As mentioned 
above, the rubber-elastic model of Ref. Q assumes that 
the layer normal and the director are rigidly locked and 
thus inevitably produces 9 = for any value of A zz - 
From their X-ray data, Nishikawa and Finkelmann con- 
clude that there is no relative tilt between layer normal 
and the director below and above the threshold. The 
Lagrangian model of Ref. Q predicts nonzero but small 
above the threshold, like in the dashed (red) curve in 
Fig. [9l Given that the angle-resolution in the experiment 
of Nishikawa and Finkelmann was of the order of one de- 
gree, it is possible that there was Sm-C like tilt in the 
experiment that was too small to detect if 6 followed a 
curve similar to the dashed (red) curve in Fig. [9] In this 
case, there is no contradiction between the experimen- 
tal data and the theoretical findings of Ref. Q and the 
present paper. 



2. Compression (A zz < 1) 

For the compression case we again use the free energy 
expression of Eq. §§§ since in the absence of a A zx there 
is no layer rotation. Compression along the layer normal 
is resisted by the layer spacing potential. There is also 
a nematic rubber elastic penalty for the chain compres- 
sion. Nematic elastomers in theory and experiment fl8j 
are known to reduce their elastic energy on compression 
along the director by rotation of the director, thereby 
presenting a shorter dimension of their chain distribution 
to the imposed strain. We find an analogous effect here. 
Numerical results for the compression case are shown in 
Fig. [H3 and illustrated in Fig. [TT1 

There is again a threshold, A c , at which rotation starts 
in the weak anchoring model. The threshold obeys the 
following polynomial: 



cr A c + 6rA c 2 (A c - 1) + (r - 1)(A' 3 - 1) = 0. 



(38) 



There is always a solution to this equation, that is, in 
principle an instability against compression should al- 
ways exist. For B 3> a t , the solution is 



A r 



1-at/B, 



(39) 



which is identical to the locked layer normal case. 



which provides an important test for the magnitude of 
the parameter at- For example if at ~ [i and B ~ 60/i 
then A c ~ 0.983, i.e. a compression of a few percent. 
For large c the instability moves to (unphysically) large 
compressions, A c ~ 0. At small deformations, before any 
rotational instability, the behavior is exactly as in the 
elongation case, given in Eq. (|3"T|) . 
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FIG. 10: The nominal shear stress, rotation angle, and defor- 
mation tensor components for an imposed Azz deformation. 
The key is as in Fig. [2] Note that the stress-deformation 
curve features a slight negative slope that is hard to see. 
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FIG. 11: An illustration of the rotation of the director, and 
the sympathetic shears for an imposed A zz < 1 compression. 

The stress-deformation curve shown in Fig. [10] is again 
non-monotonic, so is unstable in the region of negative 
slope. This unphysical feature can be resolved by re- 
sorting to a Maxwell construction as explained in Sec. 
IIIIBI As in Sec. IIIIB| this construction leads to the pre- 
diction of a plateau in the stress-deformation curve and 
microstructure that features, unless boundary conditions 
prevent this, a mixture of sheared and compressed re- 
gions. 

Reference Q focused on stretching along the layer nor- 
mal as in the experiments of Nishikawa and Finkclmann, 
and its authors chose not to consider compression along 
the layer normal. However, from the equations of Sec. Ill 
of Ref. , it is straightforward to see that the Lagrangian 
model also predicts an instability towards shear for com- 
pression along the layer, u zz < 0, where y. = ^(A T A— £) is 

the usual strain tensor. The thresholds u zz for the onset 
of u xz shear are determined by the values of u zz for which 
the effective modulus rn(u zz ) in the tilt energy density 
passes through zero from positive to negative. For the 
precise definition of the effective modulus, which is re- 
lated to the parameter at defined via Eq. ([5]) , we refer to 
Eq. (3.10b) of Ref. 0. The equation r R (u zz ) = has 2 
solutions, the one given in Eq. (3.11) of Ref. 0, and a 
corresponding one with the minus in front of the square 
root replaced by a plus. The solution with the plus per- 
tains to compression and was, therefore, not discussed in 
Ref. ML Based on the available experimental data for 
fluid [19| and elastomeric smectics [2fJ, the Lagrangian 



theory of Ref. produces the estimates u zz ~ —0.025 
and u zz « —0.0025, respectively, which are consistent 
with estimates for A c based Eq. (|39|) and the typical val- 
ues for a t and B quoted in Sec. [TTJ 

IV. DISCUSSION AND CONCLUSION 

We have explored the consequences of having a weakly 
anchored director in a microscopic model of a Sm-^4 elas- 
tomer for three deformations: in-plane shear, in-plane 
elongation, and deformation parallel to the layer normal. 
For the in-plane shear, it is found that the director rotates 
toward the extension diagonal and that initially the rota- 
tion angle is proportional to the amplitude of the shear 
applied. For in-plane elongation, and for compression 
and elongation parallel to the layer normal, the possibil- 
ity of rotation of the director leads to the prediction of 
instabilities of the system to director rotation. 

The stress-deformation curves predicted by our model 
for in-plane elongation and for compression parallel to the 
layer normal are non-monotonic. Since a region of neg- 
ative slope in the stress-deformation curve is unstable, 
this behavior is unphysical and will not be seen experi- 
mentally. Instead the elastomer should form a mixture of 
two different deformations, and exhibit a plateau in the 
stress-deformation curve, as illustrated in ijlll Bl In the 
case of imposed in-plane elongation disproportionation 
has not been reported experimentally; it may, however, 
be possible to engineer elastomers in which this behavior 
could be observed. 

The instability toward the development of shear for 
elongation along the layer normal has been discussed the- 
oretically in earlier papers 0,0- Here, we also predict an 
instability to rotation of the director under compression 
along the layer normal. This could be analyzed exper- 
imentally using the samples in which the director and 
layer normal are parallel to the normal to the film as 
reported in Q . The thresholds at which the different in- 
stabilities occur provide a useful way to determine the 
model parameters experimentally, and to find the value 
of the parameter at ■ 
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APPENDIX: GENERAL 2D MODEL 
FORMULATION 

We give here a practical method of calculating the free 
energy in the complex situations of extensions, shears 



and rotation. Simplification is possible because the de- 
formations occur in the xz subspace of A and are denoted 

by 



G 



A xx A xz 



(A.l) 



with a single y deformation (A yy = det[GJ 1 ) that acts 

to preserve volume. Using this notation, and assuming 
that the director, n remains in the xz plane, then the 
ncmatic free energy density is 



fel = \v { 



Tr 



T o-l 



G-io-C -L 



+ (dct£)- 2 }. (A.2) 



The layer normal can be calculated as follows 

(fi-x)xy 

kocA" T -z = (A-x)(xA-y) = — — ,(A.3) 



detG 



lA" 



lfi-x| 
detG 



(A.4) 



director. 



n • k = 



((G-x)xy)-n . /(n x {G ■ x)) • (n x (£ • x)) 



G • x x y 



1 - n ■ Gx 



|G-x| 



1/2 



(A.5) 



where Gx denotes the unit vector G_x/|G_x|. Combining 

the above results produces the following total free energy 
density expression 



T o-l 



G-£o-G ■£ 



(detG) 5 



\B (detG/ |G-x| 2 - (n-G-x) 



1/2 



ia t (n ■ Qjx) 



(A.6) 



The last term we require is the dot product of n and 
k to calculate the angle between the layer normal and 



This expression no longer involves the cofactor of the 
deformation gradient. 
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